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I.  INTRODUCTION 


This  report  describes  a  numerical  experiment  conducted  to  investigate  the  constraints 
placed  on  smoke  modeling  by  a  common  assumption.  That  assumption,  known  as  separability 
[  1 1.  ignores  the  differential  diffusion  properties  of  a  polydisperse  aerosol  for  simplicity. 

In  this  experiment,  the  consequences  of  this  assumption  were  examined.  Idealized, 
spherical  aerosols  of  known  source  particle  size  distribution  were  allowed  to  differentially 
diffuse  in  a  Fickian  manner.  To  avoid  geometric  confusion,  spherical  symmetry  was  assumed. 
Additionally,  the  separability  assumption  was  applied  for  average  particle  sizes  corresponding 
to  the  first  three  moments  of  the  source  particle  size  distribution.  Calculations  were  then  made 
at  selected  points  and  times  within  the  aerosol  cloud  of  total  concentration,  extinction 
coefficient,  optical  depth,  and  transmission. 

The  results  of  these  calculations  were  then  compared  to  assess  the  impact  of  the 
separability  assumption  on  the  accuracy  of  smoke  modeling.  The  theory  of  the  experiment  is 
developed  in  Sections  II  and  III.  The  conduct  of  the  experiment  is  described  and  the  results 
presented  and  discussed  in  Section  IV.  Conclusions  are  presented  in  Section  V. 

This  experiment  is  part  of  an  on-going  examination  of  the  effect  of  battlefield 
obscuration  on  the  performance  of  missile  weapons  systems  as  part  of  the  Concepts  Analysis 
and  Validation  work  area  of  the  A2I4  Missile  Technology  Program.  The  results  of  this 
experiment  will  be  used  in  the  formulation,  analysis,  and  evaluation  of  present  and  conceptual 
missile  weapon  systems. 

II.  TRANSPORT  AND  DIFFUSION 

The  transport  and  diffusion  of  a  smoke  cloud  is  usually  modeled  under  the  assumption 
of  separability  1 1 1.  This  assumption  allows  the  particle  concentration  —  size  distribution 
function  N(r,R)  where  r  is  the  particle  radius  and  R  is  position  and  time  (four  vector)  to  be 
separated  in  the  form 

N(r,R)  =  n  (r)  C(R)  (1) 

where:  n(r)  =  particle  size  distribution  function,  and 

C'(R)  =  particle  concentration  function. 


Normally.  N(r.R  i  has  units  of  number  pei  volume  per  radius  while  n(r)  has  units  of  number 
per  mass  per  radius,  and  C'(R)  has  units  of  mass  per  volume. 


The  extinction  coefficient  a  (R)  (units  of  reciprocal  length)  has  the  form 


a  (R)  = 


CO 

/N 


(r,R)a(r)  dr. 


(2) 


where  a(r)  is  the  extinction  cross  section,  usually  calculated  using  Mie  theory  [2],  In  practice, 
E  quation  (2)  is  approximated  as 


a  (R)  =  J  N(r,R)a(r)  dr 


(3) 


since  only  those  particles  of  appreciable  extinction  cross  section  or  hydrodynamic  size  are 
measurable. 


Equation  ( 3)  may  be  combined  with  Equation  ( 1 )  to  yield  the  average  <a(R)>, 


<a(R)>  =  /  n  (r)  C(R)o(r)  dr 

r. 


(4a) 


=  C  ( R)  a 


(4b) 


where: 


/ 


n(r)a(r)  dr. 


This  is  the  profit  of  the  separability  approximation;  the  specific  cross  section,  a,  is  a  constant 
within  the  cloud. 


The  particle  concentration  function  C(  R)  is  frequently  modeled  by  trivariate  gaussian 
function  that  are  approximate  solutions  of  the  diffusion  equation, 


dC 

dt 


dv2c 


(5) 
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Equation  (5)  is  appropriate  for  isotropic  diffusion  along  streamlines  and  is  much  simpler  than 
the  diffusion  equation  normally  used.  For  the  purposes  of  this  experiment,  however,  it  will 
suffice.  The  parameter  D,  the  diffusion  coefficient,  is  a  function  of  several  factors,  including 
the  particle  radius.  (See  Reference  3  )  The  form  of  D  is 


D  = 


kTg 

6irur 


(6) 


where 


k  =  Boltzman's  constant, 

T  =  absolute  temperature, 
p  =  viscosity  of  air,  and  » 

g  =  1  +  l/x  ( A^  +  A2  exp (-2A3r/£) ) 
i  =  mean  free  path  of  air  molecules 
A^  =  empirical  constants 
A1  =  1,257, 

A2  =  0.400,  and 
A3  =  0.500. 


(7) 


Equation  (7)  may  be  seen  to  be  a  function  of  temperature  and  particle  radius.  Temperature  is 
commonly  considered  in  the  modeling  of  smoke;  particle  radius  is  not.  It  is  the  purpose  of  this 
experiment  to  investigate  the  effect  of  this  neglect  on  the  properties  of  an  aerosol  cloud. 

It  may  be  shown  f  4]  that  the  form  of  the  particle  concentration  -  size  distribution  that 
satisfies  Equation  (5)  is 


N(r,R) 


G(r,R-Ro)  S(r,RQ) 


(8) 


where:  S(r,RQ)  =  source  of  particles, 

G(r,R)  =  Green  function  satisfying  Equation  (5),  and 

R  =  coordinate  of  source, 

~o 

if  no  particle  growth  occurs.  If  I  R-R<,  I  »l  R..I.  S(r.R„)  —  n(r)  S(R„).  and  Equation  (8) 
becomes 


N(r,R)  -  n  (r )  G  (r ,  R)  . 


(9) 
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The  separable  approximation  may  then  be  made  if  there  exists  an  average  radius  <r>,  such 
that 


|G(<r>  , R)  -  /  dr  Gl^Rj/fr^  -  r2)|<<  e  (10) 

rl 

where  *  is  an  experimental  measurement  error,  at  arbitrarily  many  R.  This  is  true  when 

a?  <<:  ? 


HD 


w  here  r  is  a  characteristic  time  analogous  to  the  collision  time  in  classical  transport  theory  |4). 
The  condition  expressed  by  Equation  (11)  is  commonly  assumed. 


In  this  experiment.  Equation  (9)  is  assumed:  that  is,  a  point  source  of  particles  is 
present.  Additionally,  the  medium  in  which  the  diffusion  occurs  is  assumed  to  be  isotropic  and 
steady  state.  All  coordinates  will  be  taken  relative  to  the  source  streamlines.  For  convenience, 
we  shall  assume  <r>  is  given  by 


<r> 


r2 

=  £  y*  rm  n  (r)  dr  J 


1 

m 


(12) 


where  m  will  vary  from  I  through  3.  The  average  radius  is  thus  a  moment  of  the  source  particle 
size  distribution  function. 

One  of  the  parameters  that  will  be  examined  in  this  experiment  is  the  volume  ot 
particles.  The  volume  of  particles  per  unit  volume  for  the  differential  diffusion  may  be 
calculated  from  Equation  (9)  as 


V  (R)  =  4tt/3  /  n  (r)  G  (r ,  R)  r  dr. 


(13) 


The  volumes  corresponding  to  average  radii,  called  herein  the  average  volume,  are 


<V>  =  4tt/3  G  (<r>  ,R)  f  n(r)  dr. 


(14) 
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I  he  calculation  ol  liquation  I  1.1)  will  be  further  discussed  in  the  next  section. 


Ill  OPTICAL  PROPERTIES 

In  the  previous  section,  the  extinction  coefficient  was  given  as 

C2 

a(R)  -  J  N(r,R)o(r)  dr  (3) 

rl 

which  may  be  rewritten  in  the  form 

r2 

a  (R)  -  f  n(r)  G(r,R)o(r)  dr  (15) 

rl 

by  substitution  of  Equation  (9).  We  now  introduce  the  approximation  that 

|  G  (r  +  6r,R)  -  C(r,R)|<<e  (16) 


and  that 

|o(r  +  <5r)  -  a(r)|<<iT  r2,  (17) 

for  suitably  small  6r.  This  approximation  allows  a(R)  to  be  further  approximated  as 

r  +  6r 

N  n  2 

a  (R)  -  l  ft  Q(r  )  G(r  ,R)  /  r  n(r)  dr  (18) 

~  n  n  - 


where:  6r  =  < r-  -  ri )  (N-l),  and 

Q  =  Mie  scattering  efficiency  factor". 

Another  parameter  that  will  be  examined  is  optical  depth,  the  integral  of  the  extinction 
coefficient  along  a  linc-of-sight  I  he  total  optical  depth  can  be  calculated  from  Equation  ( 19)  as 
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(19) 


( +u>  \ 

(  rn  t  6r  ) 

T  -  u  E  Q(rn)  |  f  G(rn*R)dxj 

|  J  r2n  (r)  dr .j 

The  last  parameter  to  be  examined  is  the  transmission,  which  is  given  by  the  Beer-Lambert- 
Bougert  l  aw  if  scattering  is  neglected, 

T  -  exp  (-  t)  .  (20) 

Within  this  same  approximation,  other  previously  defined  quantities  may  be 
approximated.  The  volume  of  particles.  Equation  (13).  may  be  rewritten  as 


r  + 
n 


V  (R)  -  4  it/3  l  G(rn,R) 


f  r3  n (r)  dr. 


n 


(21) 


The  average  optical  depth  and  average  transmission  corresponding  to  the  average  radii 
are 


+  0O 

< r >  -  a  /  G«r>,R)  dx,  (22) 

—  CO 


and 


<T>  =  exp  (-  <t>)  , 
respectively. 

It  should  be  noted  that  the  terms  called  averages  are  not  truly  averages  over  the  cloud. 
Rather  they  refer  only  to  the  average  radii. 

IV.  EXPERIMENTAL 

In  this  section,  the  use  of  the  equations  developed  in  Sections  11  and  111  will  be 
discussed.  Nine  excursions  will  be  described.  The  particle  size  distribution  function  n(r)  and 
the  moment  m  will  be  varied.  Three  size  distribution  functions  will  be  used,  two  Junge 
distributions  [5],  given  by 
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n (r)  =  Ar 


(23) 


-4 


and 


n<r)  =  Br"3  (24) 

where:  A  and  B  =  constants,  and  a  Deirmendjian  distribution  [2]  given  by 

n(r)  =  a  r  Y  exp(-  br)  (25) 

where,  a,  y ,  and  b  are  constants.  The  parameters  are  adjusted  to  give  the  same  total  number 
of  particles  and  for  the  Deirmendjian  distribution  to  lie  approximately  between  the  two  Junge 
distributions  for  large  r. 

Three  moments  (m  =  1 ,  2,  and  3)  were  used  corresponding  to  the  radius,  the  area,  and 
the  volume  moment.  Calculations  were  then  performed  using  these  particle  size  distributions 
of  the  volume  of  particles,  the  extinction  coefficient,  the  optical  depth,  and  the  transmission. 
The  results  of  these  calculations  are  shown  in  Figures  I  through  81. 

The  curves  of  volume  are  given  in  Figures  I  through  12.  In  all  cases,  the  volume  of 
particles  predicted  by  the  average  diffusion  parametrization  is  less  and  more  sharply 
contained  than  that  predicted  by  the  exact  formalism.  If  an  average  diffusion  parametrization 
were  used  to  analyze  experimental  data,  excessively  large  diffusion  coefficients  would  be 
required.  Exactly  similar  behavior  is  also  demonstrated  by  the  extinction  coefficient.  Figures 
13  through  24.  and  by  optical  depth.  Figures  25  through  36. 

The  surprise  in  the  experiment  may  be  found  in  the  transmission.  Figures  37  through 
81.  Large  differences  may  be  seen  in  the  transmission  curves  of  the  Junge  distributions. 
Figures  37  through  66.  although  the  clouds  are  sharp.  In  the  transmission  of  the  Deirmendjian 
haze,  however,  surprisingly  good  agreement  may  be  seen.  Further,  the  degree  of  agreement 
does  not  appear  to  vary  strongly  with  moment.  This  is  probably  due  to  the  fact  that  the 
moments  of  this  distribution  do  not  vary  rapidly  with  change  in  moment.  Additionally,  it  is 
heartening  since  this  distribution  is  most  similar  to  the  log  normal  distribution  normally  used 
in  modeling  smoke. 
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V. 


CONCLUSIONS 


I  he  volume  ol  particles.  extinction  coefficient,  optical  depth,  and  transmission  for 
aerosol  clouds  were  calculated  using  both  particle  si/e  averaged  and  si/e  dependent  diffusion. 
Non-diffusive  contributions  to  the  particle  si/e  distributions  were  not  considered,  l.arge 
differences  between  the  average  and  exact  calculations  were  evidenced  for  all  calculations  and 
particle  si/e  distributions  except  for  transmission  through  a  Deirmendjian  ha/e  distribution. 

1  he  results  of  this  experiment  do  yield  some  insight  into  the  utility  of  the  average 
diffusion  approximation  commonly  used  in  smoke  modeling.  These  models  should  be 
acceptably  accurate  for  simulation  of  transmission,  but  possibly  not  for  simulation  of 
concentration.  By  the  same  token,  these  models  should  be  validated  against  transmission 
measurements  and  not  against  concentration  measurements.  Additionally,  selection  of  an 
average  particle  si/e  distribution  is  critical. 
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Figure  1.  Average  volume  versus  distance,  first  moment  average,  n(r)  =  Mr*. 


Figure  3.  Average  volume  versus  distance,  third  moment  average,  n(r)  =  Mr*. 


Figure  5.  Average  volume  versus  distance,  first  moment  average,  n(r)  =  1/r3. 


Figure  7.  Average  volume  versus  distance,  third  moment  average,  n(r)  =  1/r3. 


Figure  10.  Average  volume  versus  distance,  second  moment  average, 
n(r)  =  Deirmendjian  haze  H. 
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Figure  11.  Average  volume  versus  distance,  third  moment  average, 
n(r)  =  Deirmendjian  haze  H. 


Figure  12.  Volume  versus  distance,  n(r)  =  Deirmendjian  haze  H. 


Figure  13.  Average  extinction  coefficient  versus  distance,  first  moment  average, 
n(r)  =  Mr*. 


Figure  1 4.  Average  extinction  coefficient  versus  distance,  second  moment  average, 
n(r)  =  1/r". 


Figure  15.  Average  extinction  coefficient  versus  distance,  third  moment  average, 
n(r)  =  1/r*. 
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Figure  16.  Extinction  coefficient  versus  distance,  n(r)  =  Mr*. 
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Figure  17.  Average  extinction  coefficient  versus  distance,  first  moment  average, 
n(r)  =  1/r3. 


Figure  1 8.  Average  extinction  coefficient  versus  distance,  second  moment  average, 
n(r)  =  1/r3. 


Figure  21.  Average  extinction  coefficient  versus  distance,  first  moment  average, 
n(r)  =  De irmendjian  haze  H. 


Figure  22.  Average  extinction  coefficient  versus  distance,  second  moment  average, 
n(r)  =  Deirmendjian  haze  H. 
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Figure  24.  Extinction  coefficient  versus  distance,  n(r)  =  Deirmendjian  haze  H. 
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Figure  25.  Average  optical  depth  versus  distance,  first  moment  average, 
n(r)  =  1/r«. 
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n(r)  =  1/r*. 


Figure  27.  Average  optical  depth  versus  distance,  third  moment  average, 
n(r)  =  1/r\ 


Figure  28.  Optical  depth  versus  distance,  n(r)  =  1/r\ 


Figure  29.  Average  optical  depth  versus  distance,  first  moment  average, 
n(r)  =  1/r3. 


Figure  30.  Average  optica)  depth  versus  distance,  second  moment  average. 
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Figure  32.  Optical  depth  versus 
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Figure  34.  Average  optical  depth  versus  distance,  second  moment  average, 
n(r)  =  Deirmendjian  haze  H. 


re  35.  Average  optical  depth  versus  distance,  third  moment  average, 
n(r)  =  Deirmendjian  haze  H. 
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Figure  39.  Transmission  versus  distance,  first  moment  average,  time  T3, 
n(r)  =  T/r*. 


Figure  41.  Transmission  versus  distance,  first  moment  average,  time  T5. 
n(r)  =  1/r*. 


Figure  42.  Transmission  versus  distance,  second  moment  average,  time  T-j, 
n{r)  =  1/r4. 
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Figure  43.  Transmission  versus  distance,  second  moment  average,  time  T2 
n(r)  =  1/r4. 


Figure  44.  Transmission  versus  distance,  second  moment  average. 
n(r)  =  1/r4. 
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Figure  45.  Transmission  versus  distance,  second  moment  average,  time  T4 
n(r)  =  1/r*. 
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igure  48.  Transmission  versus  distance,  third  moment  average,  time  T2 


Figure  53.  Transmission  versus  distance,  first  moment  average,  time  T2 
n(r)  =  1/rJ. 
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Figure  55.  Transmission  versus  distance,  first  moment  average,  time  T4, 
n(r)  =  1/r’. 


Figure  56.  Transmission  versus  distance,  first  moment  average,  time  T5, 
n(r)  =  1/r3. 


Figure  57.  Transmission  versus  distance,  second  moment  average,  time  Ti 
n(r)  =  1/r3. 


Figure  59.  Transmission  versus  distance,  second  moment  average,  time  T3 
n(r)  =  1/r3. 
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Figure  60.  Transmission  versus  distance,  second  moment  average,  time  T4, 
n(r)  =  1/r3. 
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Figure  63.  Transmission  versus  distance,  third  moment  average,  time  T2 
n(r)  “  1/r3. 


n(r)  =  1/r3. 
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Figure  65.  Transmission  versus  distance,  third  moment  average,  time  T4, 
n(r)  =  1/r3. 
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Figure  67.  Transmission  versus  distance,  first  moment  average,  time  T-j 
n(r)  =  Defrmendjian  haze  H. 
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Figure  68.  Transmission  versus  distance,  first  moment  average,  time  T2, 
n(r)  =  Deirmendjian  haze  H. 
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Figure  69.  Transmission  versus  distance,  first  moment  average,  time  13, 
n(r)  =  Deirmendjian  haze  H. 


71.  Transmission  versus  distance,  first  moment  average,  time  T5, 
n(r)  =  Deirmendjian  haze  H. 
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Figure  73.  Transmission  versus  distance, 
n(r)  =  Deirmendji an  haze  H. 
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